The evolutionary maintenance of ancient recombining sex chromosomes in the ostrich.

Sex chromosomes have evolved repeatedly across the tree of life and often exhibit extreme size dimorphism due to genetic degeneration of the sex-limited chromosome (e.g. the W chromosome of some birds and Y chromosome of mammals). However, in some lineages, ancient sex-limited chromosomes have escaped degeneration. Here, we study the evolutionary maintenance of sex chromosomes in the ostrich (Struthio camelus), where the W remains 65% the size of the Z chromosome, despite being more than 100 million years old. Using genome-wide resequencing data, we show that the population scaled recombination rate of the pseudoautosomal region (PAR) is higher than similar sized autosomes and is correlated with pedigree-based recombination rate in the heterogametic females, but not homogametic males. Genetic variation within the sex-linked region (SLR) (π = 0.001) was significantly lower than in the PAR, consistent with recombination cessation. Conversely, genetic variation across the PAR (π = 0.0016) was similar to that of autosomes and dependent on local recombination rates, GC content and to a lesser extent, gene density. In particular, the region close to the SLR was as genetically diverse as autosomes, likely due to high recombination rates around the PAR boundary restricting genetic linkage with the SLR to only ~50Kb. The potential for alleles with antagonistic fitness effects in males and females to drive chromosome degeneration is therefore limited. While some regions of the PAR had divergent male-female allele frequencies, suggestive of sexually antagonistic alleles, coalescent simulations showed this was broadly consistent with neutral genetic processes. Our results indicate that the degeneration of the large and ancient sex chromosomes of the ostrich may have been slowed by high recombination in the female PAR, reducing the scope for the accumulation of sexually antagonistic variation to generate selection for recombination cessation.

Effect of Different Types of Sequence Data on Palaeognath Phylogeny.

Palaeognathae consists of five groups of extant species: flighted tinamous (1) and four flightless groups: kiwi (2), cassowaries and emu (3), rheas (4), and ostriches (5). Molecular studies supported the groupings of extinct moas with tinamous and elephant birds with kiwi as well as ostriches as the group that diverged first among the five groups. However, phylogenetic relationships among the five groups are still controversial. Previous studies showed extensive heterogeneity in estimated gene tree topologies from conserved nonexonic elements, introns, and ultraconserved elements. Using the noncoding loci together with protein-coding loci, this study investigated the factors that affected gene tree estimation error and the relationships among the five groups. Using closely related ostrich rather than distantly related chicken as the outgroup, concatenated and gene tree-based approaches supported rheas as the group that diverged first among groups (1)-(4). Whereas gene tree estimation error increased using loci with low sequence divergence and short length, topological bias in estimated trees occurred using loci with high sequence divergence and/or nucleotide composition bias and heterogeneity, which more occurred in trees estimated from coding loci than noncoding loci. Regarding the relationships of (1)-(4), the site patterns by parsimony criterion appeared less susceptible to the bias than tree construction assuming stationary time-homogeneous model and suggested the clustering of kiwi and cassowaries and emu the most likely with ∼40% support rather than the clustering of kiwi and rheas and that of kiwi and tinamous with 30% support each.

Analysis of Sex Chromosome Evolution in the Clade Palaeognathae from Phased Genome Assembly.

Birds in the clade Palaeognathae, excluding Tinamiformes, have morphologically conserved karyotypes and less differentiated ZW sex chromosomes compared with those of other birds. In particular, the sex chromosomes of the ostrich and emu have exceptionally large recombining pseudoautosomal regions (PARs), whereas non-PARs are classified into two strata according to the date of their origins: stratum 0 and stratum 1 (S1). However, the construction and analysis of the genome sequences in these regions in the clade Palaeognathae can be challenging because assembling the S1 region is difficult owing to low sequence diversity between gametologs (Z-linked and W-linked sequences). We addressed this issue by applying the Platanus-allee assembler and successfully constructed the haplotype-resolved (phased) assembly for female emu, cassowary, and ostrich using only sequence read data derived from the Illumina platform. Comparative genomic and phylogenetic analyses based on assembled Z-linked and W-linked sequences confirmed that the S1 region of emu and cassowary formed in their common ancestor. Moreover, the interspersed repetitive sequence landscapes in the S1 regions of female emu showed an expansion of younger repetitive elements in the W-linked S1 region, suggesting an interruption in homologous recombination in the S1 region. These results provide novel insights into the trajectory of sex chromosome evolution in the clade Palaeognathae and suggest that the Illumina-based phased assembly method is an effective approach for elucidating the evolutionary process underlying the transition from homomorphic to differentiated sex chromosomes.

Chromosome-level assembly reveals extensive rearrangement in saker falcon and budgerigar, but not ostrich, genomes.

BACKGROUND: The number of de novo genome sequence assemblies is increasing exponentially; however, relatively few contain one scaffold/contig per chromosome. Such assemblies are essential for studies of genotype-to-phenotype association, gross genomic evolution, and speciation. Inter-species differences can arise from chromosomal changes fixed during evolution, and we previously hypothesized that a higher fraction of elements under negative selection contributed to avian-specific phenotypes and avian genome organization stability. The objective of this study is to generate chromosome-level assemblies of three avian species (saker falcon, budgerigar, and ostrich) previously reported as karyotypically rearranged compared to most birds. We also test the hypothesis that the density of conserved non-coding elements is associated with the positions of evolutionary breakpoint regions. RESULTS: We used reference-assisted chromosome assembly, PCR, and lab-based molecular approaches, to generate chromosome-level assemblies of the three species. We mapped inter- and intrachromosomal changes from the avian ancestor, finding no interchromosomal rearrangements in the ostrich genome, despite it being previously described as chromosomally rearranged. We found that the average density of conserved non-coding elements in evolutionary breakpoint regions is significantly reduced. Fission evolutionary breakpoint regions have the lowest conserved non-coding element density, and intrachromomosomal evolutionary breakpoint regions have the highest. CONCLUSIONS: The tools used here can generate inexpensive, efficient chromosome-level assemblies, with > 80% assigned to chromosomes, which is comparable to genomes assembled using high-density physical or genetic mapping. Moreover, conserved non-coding elements are important factors in defining where rearrangements, especially interchromosomal, are fixed during evolution without deleterious effects.

A Genetic Map of Ostrich Z Chromosome and the Role of Inversions in Avian Sex Chromosome Evolution.

Recombination arrest is a necessary step for the evolution of distinct sex chromosomes. Structural changes, such as inversions, may represent the mechanistic basis for recombination suppression and comparisons of the structural organization of chromosomes as given by chromosome-level assemblies offer the possibility to infer inversions across species at some detail. In birds, deduction of the process of sex chromosome evolution has been hampered by the lack of a validated chromosome-level assembly from a representative of one of the two basal clades of modern birds, Paleognathae. We therefore developed a high-density genetic linkage map of the ostrich Z chromosome and used this to correct an existing assembly, including correction of a large chimeric superscaffold and the order and orientation of other superscaffolds. We identified the pseudoautosomal region as a 52 Mb segment (≈60% of the Z chromosome) where recombination occurred in both sexes. By comparing the order and location of genes on the ostrich Z chromosome with that of six bird species from the other major clade of birds (Neognathae), and of reptilian outgroup species, 25 Z-linked inversions were inferred in the avian lineages. We defined Z chromosome organization in an early avian ancestor and identified inversions spanning the candidate sex-determining DMRT1 gene in this ancestor, which could potentially have triggered the onset of avian sex chromosome evolution. We conclude that avian sex chromosome evolution has been characterized by a complex process of probably both Z-linked and W-linked inversions (and/or other processes). This study illustrates the need for validated chromosome-level assemblies for inference of genome evolution.

Highly deleterious variations in COX1, CYTB, SCG5, FK2, PRL and PGF genes are the potential adaptation of the immigrated African ostrich population.

Because of variable inconvenient living conditions in some places around the world, it is difficult to collect reliable physiological data for ostriches. Therefore, this study aims to provide a comprehensive in silico insight for the nature of polymorphism of important genetic loci that are related to physiological and reproductive traits. Sixty-nine mature ostriches ranging over half of Iraq were screened. Six exonic genetic loci, including cytochrome c oxidase I (COX1), cytochrome b (CYTB), secretogranin V (SCG5), feather keratin 2-like (FK2), prolactin (PRL) and placenta growth factor (PGF) were genotyped by PCR-single stranded conformation polymorphism (SSCP). Thirty-six novel SNPs, including seventeen nonsynonymous (ns) SNPs, were observed. Several computational software programs were utilized to assess the extent of the nsSNPs on their corresponding proteins structure, function and stability. The results showed several deleterious functional and stability changes in almost all the proteins studied. The total severity of each missense mutation was evaluated and compared with other nsSNPs accumulatively. It is evident from the extensive cumulative in silico computation that both p.E34D and p.E60K in PGF have the highest deleterious effect. The cumulative predictions from the present study are an impressive guide for the genotypes of African ostriches, which bypassed the expensive protocols for wet laboratory screening, to identify the effects of variants. To the best of our knowledge, this is the first investigation of its kind on the analyses and prediction outcome of missense mutations in African ostrich populations. The highly deleterious nsSNPs in the placenta growth factor are possible adaptive mutations which might be associated with adaptation in extreme and new environments. The flow and protocol of the computational predictions can be extended for various wild animals to identify the molecular nature of adaptations.

Ancient DNA Reveals Late Pleistocene Existence of Ostriches in Indian Sub-Continent.

Ancient DNA (aDNA) analysis of extinct ratite species is of considerable interest as it provides important insights into their origin, evolution, paleogeographical distribution and vicariant speciation in congruence with continental drift theory. In this study, DNA hotspots were detected in fossilized eggshell fragments of ratites (dated ≥25000 years B.P. by radiocarbon dating) using confocal laser scanning microscopy (CLSM). DNA was isolated from five eggshell fragments and a 43 base pair (bp) sequence of a 16S rRNA mitochondrial-conserved region was successfully amplified and sequenced from one of the samples. Phylogenetic analysis of the DNA sequence revealed a 92% identity of the fossil eggshells to Struthio camelus and their position basal to other palaeognaths, consistent with the vicariant speciation model. Our study provides the first molecular evidence for the presence of ostriches in India, complementing the continental drift theory of biogeographical movement of ostriches in India, and opening up a new window into the evolutionary history of ratites.

Multiple IgH Isotypes Including IgD, Subclasses of IgM, and IgY Are Expressed in the Common Ancestors of Modern Birds.

Although evolutionarily just as ancient as IgM, it has been thought for many years that IgD is not present in birds. Based on the recently sequenced genomes of 48 bird species as well as high-throughput transcriptome sequencing of immune-related tissues, we demonstrate in this work that the ostrich (Struthio camelus) possesses a functional δ gene that encodes a membrane-bound IgD H chain with seven CH domains. Furthermore, δ sequences were clearly identified in many other bird species, demonstrating that the δ gene is widely distributed among birds and is only absent in certain bird species. We also show that the ostrich possesses two µ genes (µ1, µ2) and two υ genes (υ1, υ2), in addition to the δ and α genes. Phylogenetic analyses suggest that subclass diversification of both the µ and υ genes occurred during the early stages of bird evolution, after their divergence from nonavian reptiles. Although the positions of the two υ genes are unknown, physical mapping showed that the remaining genes are organized in the order µ1-δ-α-µ2, with the α gene being inverted relative to the others. Together with previous studies, our data suggest that birds and nonavian reptile species most likely shared a common ancestral IgH gene locus containing a δ gene and an inverted α gene. The δ gene was then evolutionarily lost in selected birds, whereas the α gene lost in selected nonavian reptiles. The data obtained in this study provide significant insights into the understanding of IgH gene evolution in tetrapods.

Development of qualitative and quantitative PCR analysis for meat adulteration from RNA samples.

Total RNA samples were used to establish qualitative and quantitative PCR-based methods for assessing meat adulteration. The primers were designed based on the mRNA sequences of troponin I (TnI), mitochondrial ribosomal protein (MRP) and tropomodulin genes to distinguish chicken, pork, goat, beef and ostrich. There was no cross reaction between the primers, and the detection limit of the cDNA template was 0.01 and 20 ng in simplex PCR and multiplex PCR, respectively. In the low temperature storage test, the detection limits of cDNA template with 10 and 1 ng were determined at 4 °C and -80 °C. In quantitative assay, the precision of real-time PCR analysis expressed as a coefficient of variation (CV) ranged from 0.25% to 5.24% and the trueness, expressed as an error, ranged from 0.28% to 6.98% for adulteration. Thus, herein, we provided alternative tools for the assessment of meat adulteration using mRNA-based PCR methods.

Molecular cloning, expression and bioactivity of B cell activating factor (BAFF) in African ostrich.

B cell activating factor (BAFF), which belongs to the tumor necrosis factor (TNF) family, is testified to play a critical role in B cell survival, proliferation, maturation and immunoglobulin secretion. In the present study, the cDNA of open reading frame (ORF) in African ostrich (Struthio camelus) BAFF (designated OsBAFF) was cloned by reverse transcription-PCR (RT-PCR). The OsBAFF gene encodes a 288-amino acid protein containing a predicted transmembrane domain and a putative furin protease cleavage site like BAFFs from chicken (cBAFF), quail (qBAFF), duck (dBAFF), goose (gBAFF) and dove (doBAFF). RT-PCR analysis showed that the OsBAFF gene is strongly expressed in the bursa of Fabricius, thymus, spleen, and bone marrow. The soluble OsBAFF had been cloned into pET28a. SDS-PAGE and Western blotting analysis confirmed that the soluble fusion protein His-OsBAFF was efficiently expressed in Escherichia coli Rosset (DE3). In vitro, purified OsBAFF was not only able to promote the survival of African ostrich bursal lymphocytes, but also able to co-stimulate proliferation of mouse splenic B cells. The expression of OsBAFF in lymphocyte cells was higher than the control after LPS stimulation. These findings indicated that OsBAFF plays an important role in survival and proliferation of African ostrich bursal lymphocytes, which may provide valuable information for research into the immune system of African ostrich and OsBAFF may serve as a potential immunologic factor for enhancing immunological efficacy in African ostrich and any other birds.

Evolutionary analysis of the female-specific avian W chromosome.

The typically repetitive nature of the sex-limited chromosome means that it is often excluded from or poorly covered in genome assemblies, hindering studies of evolutionary and population genomic processes in non-recombining chromosomes. Here, we present a draft assembly of the non-recombining region of the collared flycatcher W chromosome, containing 46 genes without evidence of female-specific functional differentiation. Survival of genes during W chromosome degeneration has been highly non-random and expression data suggest that this can be attributed to selection for maintaining gene dose and ancestral expression levels of essential genes. Re-sequencing of large population samples revealed dramatically reduced levels of within-species diversity and elevated rates of between-species differentiation (lineage sorting), consistent with low effective population size. Concordance between W chromosome and mitochondrial DNA phylogenetic trees demonstrates evolutionary stable matrilineal inheritance of this nuclear-cytonuclear pair of chromosomes. Our results show both commonalities and differences between W chromosome and Y chromosome evolution.

Effect of boric acid supplementation of ostrich water on the expression of Foxn1 in thymus.

Foxn1 is essential for thymus development. The relationship between boric acid and thymus development, optimal dose of boric acid in ostrich diets, and the effects of boric acid on the expression of Foxn1 were investigated in the present study. Thirty healthy ostriches were randomly divided into six groups: Group I, II, III, IV, V, VI, and supplemented with boric acid at the concentration of 0 mg/L, 40 mg/L, 80 mg/L, 160 mg/L, 320 mg/L, 640 mg/L, respectively. The histological changes in thymus were observed by HE staining, and the expression of Foxn1 analyzed by immunohistochemistry and western blot. TUNEL method was used to label the apoptotic cells. Ostrich Foxn1 was sequenced by Race method. The results were as following: Apoptosis in ostrich thymus was closely related with boric acid concentrations. Low boric acid concentration inhibited apoptosis in thymus, but high boric acid concentration promoted apoptosis. Foxn1-positive cells were mainly distributed in thymic medulla and rarely in cortex. Foxn1 is closely related to thymus growth and development. The nucleotide sequence and the encoded protein of Foxn1 were 2736 bases and 654 amino acids in length. It is highly conserved as compared with other species. These results demonstrated that the appropriate boric acid supplementation in water would produce positive effects on the growth development of ostrich thymus by promoting Foxn1 expression, especially at 80 mg/L, and the microstructure of the thymus of ostrich fed 80 mg/L boric acid was well developed. The supplementation of high dose boron (>320 mg/L) damaged the microstructure of thymus and inhibited the immune function by inhibiting Foxn1 expression, particularly at 640 mg/L. The optimal dose of boric acid supplementation in ostrich diets is 80 mg/L boric acid. The genomic full-length of African ostrich Foxn1 was cloned for the first time in the study.

Complex evolutionary trajectories of sex chromosomes across bird taxa.

Sex-specific chromosomes, like the W of most female birds and the Y of male mammals, usually have lost most genes owing to a lack of recombination. We analyze newly available genomes of 17 bird species representing the avian phylogenetic range, and find that more than half of them do not have as fully degenerated W chromosomes as that of chicken. We show that avian sex chromosomes harbor tremendous diversity among species in their composition of pseudoautosomal regions and degree of Z/W differentiation. Punctuated events of shared or lineage-specific recombination suppression have produced a gradient of "evolutionary strata" along the Z chromosome, which initiates from the putative avian sex-determining gene DMRT1 and ends at the pseudoautosomal region. W-linked genes are subject to ongoing functional decay after recombination was suppressed, and the tempo of degeneration slows down in older strata. Overall, we unveil a complex history of avian sex chromosome evolution.

Phenotypic and genetic variability for body weight of ostriches (Struthio camelus).

1. objective of the present study was to estimate heritability for hatch weight (HW), body weight at 90 (W90) and 180 (W180) d of age in Brazilian ostriches. 2. The heritability estimates were obtained through the restricted maximum likelihood method for an animal model. Least squares method was used to define possible fixed effects for consideration in the analysis model. 3. The estimates of heritability were 0.42 ± 0.05, 0.16 ± 0.04 and 0.24 ± 0.08 for HW, W90 and W180, respectively. 4. Heritability estimates showed that there are sufficient levels of additive genetic variation present in the traits studied, and these traits may respond to selection.

Molecular cloning, characterization and tissue distribution of two ostrich ß-defensins: AvBD2 and AvBD7.

Avian ß-defensins (AvBDs) are a family of small antimicrobial peptides that play important roles in the innate immunity of birds. Herein, we report on two new ostrich AvBD genes, AvBD2 and AvBD7, which were isolated from the bone marrow of ostriches (Struthio camelus). The coding regions of ostrich AvBD2 and AvBD7 comprised 195 bp and 201bp, which encoded 64 and 66 amino acids, respectively. Homology analysis showed that ostrich AvBD2 had the highest similarity (up to 86%) with the swan goose (Anser cygnoides) AvBD2, while ostrich AvBD7 shared the highest similarity (up to 81%) with chicken AvBD7. Analysis of the codon-usage bias showed that the two ostrich AvBDs had different codon-usage patterns from other AvBDs. The two synthetic AvBD peptides exhibited antibacterial activities against both Gram-positive and Gram-negative bacteria, and these activities decreased significantly in the presence of 100mM NaCl (P<0.01). Real-time reverse transcription-polymerase chain reaction analysis showed that AvBD2 and AvBD7 were widely expressed at different levels in 17 different tissues. This is the first report of the nucleotide sequences of ostrich AvBDs. Further investigations of these two AvBDs may help us to gain new insights into the immune defense system of the ostrich and to make subsequent therapeutic use of ostrich defensins.

Development of loop-mediated isothermal amplification (LAMP) assay for rapid and sensitive identification of ostrich meat.

Animal species identification is one of the primary duties of official food control. Since ostrich meat is difficult to be differentiated macroscopically from beef, therefore new analytical methods are needed. To enforce labeling regulations for the authentication of ostrich meat, it might be of importance to develop and evaluate a rapid and reliable assay. In the present study, a loop-mediated isothermal amplification (LAMP) assay based on the cytochrome b gene of the mitochondrial DNA of the species Struthio camelus was developed. The LAMP assay was used in combination with a real-time fluorometer. The developed system allowed the detection of 0.01% ostrich meat products. In parallel, a direct swab method without nucleic acid extraction using the HYPLEX LPTV buffer was also evaluated. This rapid processing method allowed detection of ostrich meat without major incubation steps. In summary, the LAMP assay had excellent sensitivity and specificity for detecting ostrich meat and could provide a sampling-to-result identification-time of 15 to 20 minutes.

Reference genes for quantitative gene expression studies in multiple avian species.

Quantitative real-time PCR (qPCR) rapidly and reliably quantifies gene expression levels across different experimental conditions. Selection of suitable reference genes is essential for meaningful normalization and thus correct interpretation of data. In recent years, an increasing number of avian species other than the chicken has been investigated molecularly, highlighting the need for an experimentally validated pan-avian primer set for reference genes. Here we report testing a set for 14 candidate reference genes (18S, ABL, GAPDH, GUSB, HMBS, HPRT, PGK1, RPL13, RPL19, RPS7, SDHA, TFRC, VIM, YWHAZ) on different tissues of the mallard (Anas platyrhynchos), domestic chicken (Gallus gallus domesticus), common crane (Grus grus), white-tailed eagle (Haliaeetus albicilla), domestic turkey (Meleagris gallopavo f. domestica), cockatiel (Nymphicus hollandicus), Humboldt penguin (Sphenicus humboldti), ostrich (Struthio camelus) and zebra finch (Taeniopygia guttata), spanning a broad range of the phylogenetic tree of birds. Primer pairs for six to 11 genes were successfully established for each of the nine species. As a proof of principle, we analyzed expression levels of 10 candidate reference genes as well as FOXP2 and the immediate early genes, EGR1 and CFOS, known to be rapidly induced by singing in the avian basal ganglia. We extracted RNA from microbiopsies of the striatal song nucleus Area X of adult male zebra finches after they had sang or remained silent. Using three different statistical algorithms, we identified five genes (18S, PGK1, RPS7, TFRC, YWHAZ) that were stably expressed within each group and also between the singing and silent conditions, establishing them as suitable reference genes. In conclusion, the newly developed pan-avian primer set allows accurate normalization and quantification of gene expression levels in multiple avian species.

Ancient DNA reveals elephant birds and kiwi are sister taxa and clarifies ratite bird evolution.

The evolution of the ratite birds has been widely attributed to vicariant speciation, driven by the Cretaceous breakup of the supercontinent Gondwana. The early isolation of Africa and Madagascar implies that the ostrich and extinct Madagascan elephant birds (Aepyornithidae) should be the oldest ratite lineages. We sequenced the mitochondrial genomes of two elephant birds and performed phylogenetic analyses, which revealed that these birds are the closest relatives of the New Zealand kiwi and are distant from the basal ratite lineage of ostriches. This unexpected result strongly contradicts continental vicariance and instead supports flighted dispersal in all major ratite lineages. We suggest that convergence toward gigantism and flightlessness was facilitated by early Tertiary expansion into the diurnal herbivory niche after the extinction of the dinosaurs.

Embryonic development in the ostrich (Struthio camelus) during the first 7 days of artificial incubation.

1. Early development of ostrich embryos was investigated in relation to time of egg collection and genotype. 2. A total of 321 ostrich eggs were collected during the 2008 and 2009 breeding seasons and the development of the embryo for up to the first 168 h of incubation was described and analysed. A sample of the incubated eggs was weighed and opened daily to investigate developmental changes. 3. In fresh eggs, the blastoderm contained a round, translucent dark area pellucida (AP) in the centre, with a surrounding thin white ring, likely to be the beginning of the area opaca (AO). Fresh eggs were considered infertile if the blastoderm was absent and instead numerous white droplets were present surrounded by vacuoles. 4. The average blastoderm area of a fresh fertile egg was 15.8 mm(2), increasing to 143.3 mm(2) after 2 d of incubation. By 72 h of incubation the area vasculosa (AV) was discernible in the posterior half of the blastoderm. 5. At 48 h of incubation the blastoderm area in eggs from the South African Black genotype (SAB) × Zimbabwean Blue genotype (ZB) crosses (104.5 ± 18.6 mm(2)) was lower than the pure SAB (141.0 ± 10.5 mm(2)), ZB (161.7 ± 13.5 mm(2)) and ZB × SAB crosses (166.1 ± 14.2 mm(2)). 6. Embryo length was 5.01 mm after 72 h of incubation and 14.5 mm after 168 h of incubation. At 168 h of incubation AV lengths for both ZB × SAB (53.2 mm) and SAB × ZB crosses (54.1 mm) were longer than in embryos from the pure breeds. 7. Results from this study can be put to practical use when determining whether eggs are infertile or fertile and also in investigating the age of early embryonic mortalities.

Old but not (so) degenerated--slow evolution of largely homomorphic sex chromosomes in ratites.

Degeneration of the nonrecombining chromosome is a common feature of sex chromosome evolution, readily evident by the presence of a pair of largely heteromorphic chromosomes, like in eutherian mammals and birds. However, in ratites (order Palaeognathae, including, e.g., ostrich), the Z and W chromosomes are similar in size and largely undifferentiated, despite avian sex chromosome evolution was initiated > 130 Ma. To better understand what may limit sex chromosome evolution, we performed ostrich transcriptome sequencing and studied genes from the nonrecombining region of the W chromosome. Fourteen gametologous gene pairs present on the W chromosome and Z chromosome were identified, with synonymous sequence divergence of 0.027-0.177. The location of these genes on the Z chromosome was consistent with a sequential increase in divergence, starting 110-157 and ending 24-30 Ma. On the basis of the occurrence of Z-linked genes hemizygous in females, we estimate that about one-third of the Z chromosome does not recombine with the W chromosome in female meiosis. Pairwise d(N)/d(S) between gametologs decreased with age, suggesting strong evolutionary constraint in old gametologs. Lineage-specific d(N)/d(S) was consistently higher in W-linked genes, in accordance with the lower efficacy of selection expected in nonrecombining chromosomes. A higher ratio of GC > AT:AT > GC substitutions in W-linked genes supports a role for GC-biased gene conversion in differentially driving base composition on the two sex chromosomes. A male-to-female (M:F) expression ratio of close to one for recombining genes and close to two for Z-linked genes lacking a W copy show that dosage compensation is essentially absent. Some gametologous genes have retained active expression of the W copy in females (giving a M:F ratio of 1 for the gametologous gene pair), whereas for others W expression has become severely reduced resulting in a M:F ratio of close to 2. These observations resemble the patterns of sex chromosome evolution seen in other avian and mammalian lineages, suggesting similar underlying evolutionary processes, although the rate of sex chromosome differentiation has been atypically low. Lack of dosage compensation may be a factor hindering sex chromosome evolution in this lineage.